Generic slow relaxation in a stochastic sandpile 
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Abstract 

Simulations of a stochastic fixed-energy sandpile in one and two 
dimensions reveal slow relaxation of the order parameter, even far 
from the critical point. The decay of the activity is best described 
by a stretched-exponential form. The persistence probability (for 
a site not to have toppled up to time t), also exhibits stretched- 
exponential relaxation. The results suggest a connection between 
sandpile models and structural glasses. 
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Sandpile models have attracted much interest in recent years, as paradigms 
of scale-invariance in the apparent absence of tuning parameters , and 

as intriguing examples of absorbing-state phase transitions PJ6||7|P,P, p!0| , p!l| . 
Most studies of sandpiles have focused on the scale-invariant stationary state 
under slow driving [0,0, or on scaling properties in the vicinity of the 
absorbing-state transition . Relatively little attention has been 

given to the dynamic properties of sandpiles away from the critical point. 

A central feature of sandpile models is the presence of a conserved field, 
the density of particles. This field couples to the activity density, which is the 
order parameter. When, as in the case of FES, the conserved field is frozen in 
the absence of activity, the critical behavior is expected to fall in a universality 
class distinct from that of directed percolation (DP) |]I3|. (In the absence of 
such a conservation law, DP is generic for absorbing-state phase transitions 
|T^,T^]. Relaxation of the order parameter to its quasi-stationary value in a 



sandpile at the critical density is best characterized as a power-law 0, but 
with certain anomalies in the one-dimensional case [|T0|. Besides having an 
effect on the critical behavior, it is reasonable to expect the conservation law 
to modify relaxation properties away from the critical point as well. 

In simple models exhibiting an absorbing-state phase transition, such as 
the contact process [|^, relaxation to the stationary state is expected to be 
monotonic and exponential, away from the critical point ||21|. The simulation 
results reported below show that in sandpiles, off-critical relaxation is consid- 
erably slower, following a stretched-exponential or, in certain cases, algebriac 
decay. Stretched-exponential functions have been reported for avalanche size 
distributions in experiments on sand and rice piles P^J23], and in granular 



avalanche models pj], but not, to my knowledge, in the context of sandpile 
relaxation dynamics. 

I study a variant of Manna's sandpile p5| , defined on a lattice of L'^ sites (in 
d dimensions), with periodic boundaries. The configuration is specified by the 
number of particles = 0, 1, 2, ... at each site; sites with Zi >2 are said to be 
active. A Markovian dynamics is defined by the toppling rate, which is unity 
for all active sites, and zero for sites with Zi < 2. When a site i topples, it sends 
two particles to adjacent sites {zi ^ Zi — 2); the particles move independently 
to randomly chosen nearest neighbors j and j' {j, j' G {i + — 1} in the one- 
dimensional case). Thus j = j' with probability l/2d. The dynamics conserves 
the number of particles, A^, which is fixed by the initial configuration. 

The system evolves via a sequential dynamics: the next site to topple is 
chosen at random from a list of active sites; the time increment associated with 
a toppling is At = l/N^, where is the number of active sites just prior 
to the event. Initial configurations are generated by distributing (L'^ parti- 
cles randomly over the lattice, yielding an initial distribution that is spatially 
homogeneous and uncorrelated. Previous studies confirm the existence of a 
continuous phase transition from an absorbing to an active phase at a particle 
density Cc = 0.94885 in one dimension [jlO[, and 0.71695 in 2-d 0. For ( < (c, 



the stationary value of the order parameter (the density of active sites) is zero. 
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I begin with the simplest case, relaxation of the order parameter p in one 
dimension, for C, < Cc- For ( = 0.5, far below the critical value, the pattern of 
relaxation is essentially the same for systems of 1000 or more sites. Figure 1 
shows that p{t) (obtained from averages over 5 x 10^ realizations of a system of 
5000 sites), is well described by a stretched exponential, p{t) ~ exp[— (t/to)^] 
with f3 = 0.45 (to represents a characteristic timescale for relaxation). In this 
and subsequent analyses, the exponent (3 is determined using the criterion of 
zero curvature, in the asymptotic region of the graph of ln[p(t) — p] versus t'^ 
{p is the asymptotic activity density, which is of course zero for ^ < ^c). A 
rough estimate of the uncertainty in j3 values is ±0.02. 

Next I examine the relaxation for ( = 0.9, much nearer the critical point. 
There are now considerable finite-size effects and much larger systems (up to 
10^ sites) are required to observe the asymptotic behavior. For small systems 
p{t) appears to decay faster than a stretched exponential, but as L is increased, 
the long-time slow relaxation develops. For L = 10^, p{t) again follows a 
stretched exponential, with /3 = 0.39. It should be emphasized that while 
the simulation data cannot be interpreted as proving stretched-exponential 
relaxation, they do serve to rule out definitively both exponential and power- 
law relaxation of p{t) in the subcritical regime. 

In the supercritical regime, the relaxation of the order parameter to its 
stationary value p is nonmonotonic. p{t) decays rapidly at first, and then slowly 
approaches the stationary value from below. For = 1, I find p = 0.118222(8), 
with Ap = p — pit) decaying asymptotically as a power-law, ~ t~'^-^'^ (see Fig. 
2). The initial decay is again well described by a stretched exponential, with 
an exponent /3 = 0.28. For somewhat higher densities {( = 1.2 and 1.3), the 
asymptotic approach to p (from below), is again via a power law, with an 
exponent of 0.5 - 0.52. 



I also studied a variant of the stochastic sandpile introduced in Ref. 
in which the toppling rate at a site with z particles is z{z — l). (In this case 
(c = 0.9493(2).) For p = 1/2, the decay to p = is described with high 
precision by a stretched exponential with (3 = 0.475(25). In the supercritical 
regime the approach to the stationary value is nonmonotonic, similar to that 
observed above. Unlike the constant-rate model, however, the relaxation is 
well described by an expression consisting of two stretched exponentials: 

pit) =p + Aiexp[-Xit'^^] - A2exp[-X2t^/^], (1) 

with Ai, A2, Ai, and A2 all positive constants (Ai ^ A2, and Ai ^ A2). 

Slow relaxation appears to be robust under changes in the toppling rate. 
Of equal interest is the nature of relaxation in two or more dimensions. At a 
density of C = 0.5 (well below Cc = 0.71695(5)), studies using lattices of 
sites, with L = 40, 80, and 160, reveal that the order parameter decays as 
a stretched exponential, but with an exponent (3 = 0.81. Nearer the critical 
point i( = 0.7), finite-size effects are prominent, as was found in the one- 
dimensional case. Studies using system sizes L of up to 1280 again confirm 
stretched-exponential decay with (3 = 0.81, as shown in Fig. 3. 
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The relaxation above the critical density (in a system of size L = 160 at C = 

0. 75; p = 0.04995), is nonmonotonic, as in one dimension. The excess density 
is again well described by the difference between two stretched exponentials, 
as in Eq. (|l]), but in this case the exponent (3 associated with the early decay 
is about 0.4, while the long-time approach to p (from below) is characterized 
by an exponent of 0.8 (see Fig. 3, inset). 

As we have seen, relaxation of the order parameter in the stochastic 
sandpile is characterized by stretched-exponential (or in some cases, power- 
law) functions. Recently, O'Donoghue and Bray demonstrated stretched- 
exponential decay of the persistence probability, i.e., that a given site has 
never been visited by a diffusing particle, in certain one-dimensional reac- 
tion/diffusion processes. This suggests a study of persistence in the sandpile. 
In fact, since a site that is initially below threshold for toppling must be vis- 
ited by another diffusing particle or particles before it can topple, persistence 
seems particularly relevant to sandpile relaxation dynamics. (Studies of dy- 
namic critical exponents for persistence in conserved lattices gases, which are 
closely related to sandpiles, were reported by Liibeck [^]; here, however, we 
focus on the dynamics away from the critical point.) 

In studying persistence in a sandpile, it appears useful to group sites ac- 
cording to their initial occupation number z. Let p{t; z) denote the persistence 
probability of sites whose initial height is z. For 2; = or 1, I define p{t\ z) as 
the probability that a site has not toppled up to time t. For z > 2, however, 
such a definition is not very interesting, since (in the constant-rate model), the 
probability not to have toppled up to time t is simply e~*. For z >2, therefore, 
p{t\ z) is defined as the probability not to have toppled a second time. I study 
the one-dimensional fixed-rate model with particle densities C, = 0.5, 0.9, and 

1, as above. 

For ( < (c, all activity ceases after a finite time, so that the persistence 
probability approaches a nonzero value p{z) at long times. For ( = 0.5, the 
asymptotic relaxation of p{t; z) to p(z) again follows a stretched exponential; 
the best estimates for the exponent [3 are 0.42 for z = and (3 = 0.39 for z = 1, 
2 and 3. For C = 0.9, a study using L = W yielded /3 ~ 0.44, 0.47, 0.49, and 
0.41 for z = 0, 1, 2, and 3, respectively. Finally, for = 1, the asymptotic 
decay of the persistence probability (studied on a ring of 2 x 10*^ sites) follows a 
stretched exponential with j3 values of 0.42, 0.41, 0.38, and 0.45 for the various 
z values. The persistence probabilitites for ( = 0.9 and C = 1 ^ire shown in Fig. 
4. (In all of these studies, the exponents for z = 2 and 3 are less certain, due 
to poorer statistics, and to the fact that stretched-exponential behavior sets 
in later than for z = or 1. I find no clear trend in the estimated exponent 
value, either with z or with (.) The foregoing results can be summarized as 
indicating that in one dimension, the asympotic relaxation of the persistence 
probability p{t; z) follows a stretched exponential with an exponent f3 of about 
0.42. 

The simulation results indicate that the relaxation of the order parameter 
and of the persistence probabihty, in a stochastic sandpile away from its crit- 
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ical point, generically exhibits stretched-exponential scaling. The exponent 
f3 generally takes values in the neighborhood of 0.4 in the one-dimensional 
case, while in two dimensions (for which only order-parameter relaxation was 
studied) the value is about 0.8. Stretched-exponential relaxation may be un- 
derstood as a consequence of the depletion of active elements, leading to a 
decreasing relaxation rate, {1/ p)\dp/dt\, as time goes on. The larger exponent 
in two dimensions is plausible, given the larger number of paths in configura- 
tion space. Detailed explanations of the streteched-exponential form, and of 
specific exponent values, remain as open challenges. 

The observation of stretched-exponential relaxation suggests a connection 
with glassy dynamics. Superficially, the sandpile model seems to have little 
connection with dynamics of a dense fluid, but the two problems are related 
at a more abstract level. In dense fluids, not all the empty space in the system 
is available for particle movement, leading to a highly cooperative dynamics 
p8| , |29[| . Some larger grouping of "voids" is required for relaxation on scales 
beyond that of the local cage. (For simplicity, I frame the analogy in terms 
of a hard-sphere fluid, in which the relevant variable is density not temper- 
ature.) Similarly, in a sandpile not all particles are available for movement 
(and relaxation): only those with companions in the same cell are mobile. 
These observations suggest a parallel between particles in the sandpile and 
parcels of unoccupied space in a dense fluid. The dynamic arrest in the fluid, 
when the density approaches an (apparent) transition value, corresponds to 
the absorbing-state phase transition in the sandpile. The associated order pa- 
rameter - mobile free volume in the fluid - vanishes in the absorbing phase. 
(The total unoccupied volume is of course conserved, but the fraction that is 
mobile or "active" is not.) From this vantage, the liquid-glass transition cor- 
responds to a dynamic transition, akin to directed percolation, of free volume, 
as opposed to the static percolation transition suggested some time ago by 



Cohen and Grest |30|. While the sandpile models studied to date seem too 
simplistic to capture the dynamics of a dense fluid near the glassy state, the 
results of the present study suggest that it would be worthwhile pursuing this 
analogy by developing somewhat more realistic particle models, and studying 
diffusion and the response to external driving, in the hope of gaining further 
insight into real glasses. 
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FIGURE CAPTIONS 



FIG. 1. Relaxation of the order parameter in the one- dimensional sandpile, 
C = 0.5, L = 5000. Inset: semi-log plot of the same data. 

FIG. 2. Relaxation of the order parameter in the one- dimensional sandpile, 
C= 1, L = 2 X 10*^. Inset: log-log plot of the excess Ap = p{t) — p. 

FIG. 3. Relaxation of the order parameter in the two-dimensional sandpile, 
C = 0.7, L = 1280. Inset: ln|Ap| versus Int for ( = 0.75; points: simulation; 
solid curve: fit using two stretched exponentials as described in text. 

FIG. 4. Persistence probabilities pit; z) for 2; = 0, 1, 2, and 3 (upper to lower) 
in the one-dimensional model with C, — 1. Inset: a similar plot for C, — 0.9. 
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